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First-principles dynamical CPA (Coherent-Potential Approximation) for electron correlations has 
' been developed further by taking into account higher-order dynamical corrections with use of the 

asymptotic approximation. The theory is applied to the investigations of a systematic change of 
f— I , excitation spectra in 3d transition metals from Sc to Cu at finite temperatures. It is shown that the 

r-{ ■ dynamical effects damp main peaks in the densities of states (DOS) obtained by the local density 

approximation to the density functional theory, reduce the band broadening due to thermal spin 
fluctuations, create the Mott-Hubbard type bands in the case of fee Mn and fee Fe, and create a small 
hump corresponding to the '6 eV satellite in the case of Co, Ni, and Cu. Calculated DOS explain 
the X-ray photoelectron spectroscopy data as well as the bremsstrahlung isochromat spectroscopy 
data. Moreover, it is found that screening effects on the exchange energy parameters are significant 
^ ■ for understanding the spectra in magnetic transition metals. 

^ ' PACS numbers: 71.20.Be,78.20.Bh,75.10.Lp,78.70.En 
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I. INTRODUCTION 



, Electron correlations play an important role for understanding the electronic structure, magnetism, metal-insulator 
(~| ■ transition , and the high-temperature superconductivity in the solid-state physics, and thus a large number of theories 
O \ have been proposed to describe correlated electron systems^. Especially, in the case of magnetism, theories of electron 
CJ . correlations have been developed over fifty years to explain the ferromagnetism of transition metals, since the Hartree- 
Fock approximation was recognized to overestimate their magnetic ordering energy. 

Gutzwiller=i"i. proposed a variational theory which takes into account on-site correlations by controlling the probabil- 
ity amplitudes of doubly occupied states, and showed that electron correlations much instabilize the ferromagnetism. 
HubbardS'S developed a Green function method making use of the equation of motion method and a decoupling ap- 
proximation that leads to an alloy-analogy picture. He succeeded in describing the metal-insulator transition as well as 
, the instability of ferromagnetism due to electron correlations. Kanamorii took into account the multiple scattering of 
. • electrons in the low density limit, and showed that the effective Coulomb interaction for the ferromagnetic instability 
is extremely renormalized by electron-electron interactions. 

Above mentioned theories are limited to the ground state. Cyrot^ extended to finite temperatures an idea of 
alloy-analogy approximation for electron correlations proposed by Hubbard, on the basis of the functional inte- 
gral method^i"— . He explained the T-P phase diagram for metal-insulator transitions qualitatively. Hubbard^'^ and 
HasegawaM independently developed the single-site spin fluctuation theory using the coherent potential approximation 
(CPA) ^^i^^ . They showed that thermal spin fluctuations much reduce the Curie temperatures obtained by the Stoner 
; ^ ' theory for band calculations based on the local density approximation (LDA) to the density functional theorj*ii. 
Ci ' The single-site spin fluctuation theory reduces to the Hartree-Fock one at zero temperature because it is based 
on a high-temperature approximation, i.e., the static approximation to the functional integral method. Therefore 
the theory does not take into account the ground-state electron correlations as found by Gutzwiller, Hubbard, and 
Kanamori. Kakehashi and Fuldei^ proposed a variational theory which adiabatically takes into account such corre- 
lations at finite temperatures, and found further reduction of Curie temperature. Finally, Kakehashi proposed the 
dynamical CPA which completely takes into account the dynamical charge and spin fluctuations within the single-site 
approximation, and clarified the dynamical effects on the momentum distribution, magnetic moment as well as excita- 
tion spectra using the Monte-Carlo technique. In the next paper— which we refer to I in the following, we developed 
an analytic method to the dynamical CPA, using the harmonic approximation. In the recent paper^^ which we refer 
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to II, we proposed the first-principles dynamical CPA which combines the dynamical CPA with the tight-binding 
linear-mufhntin-orbital (TB-LMTO)^^ base LDA-I-U Hamiltonian2^. Within the 2nd-order dynamical corrections to 
the static approximation, we have shown that the dynamical CPA can describe the finite temperature properties of 
excitations and magnetism in Fe and Ni quantitatively or semiquantitatively. 

In this paper, we develop further the first-principles dynamical CPA by taking into account higher-order dynamical 
corrections within an asymptotic approximation, and investigate a systematic change of excitation spectra in 3d 
transition metals from Sc to Cu. We will clarify the dynamical effects on the excitation spectra in 3d series at finite 
temperatures, and will explain systematic change of X-ray photoelectron spectroscopy (XPS) datsM as well as the 
bremsstrahlung isochromat (BIS) data2^. 

Similar calculations have recently been performed at the ground state by Belashchenko et. al^ on the basis of 
the self-consistent local GW approximation. Their results of the first-principles calculations, however, do not well 
describe the main peak positions in the XPS and the BIS data, and seem to require further development of the theory. 

As we have proven in the separate papers^^'^^ , the dynamical CPA is equivalent to the many-body CPA^^ developed 
in the disordered system, the dynamical mean field theory (DMFT) in the metal-insulator transition in infinite 
dimensions^"—, and the projection operator method CPA (PM-CPA)^'* for excitation problems. The dynamical CPA 
was originally developed to describe the finite-temperature magnetism in metallic systems starting from the static 
approximation exact in the high-temperature limit. The theory can treat the transverse spin fluctuations for arbitrary 
d electron number at finite temperatures, though it is not easy in the traditional quantum Monte-Carlo approach 
(QMC)-- . Moreover the theory allows us to calculate excitation spectra up to the temperatures much lower than 
those calculated by the QMC because the dynamical CPA is an analytic theory which does not rely on the statistical 
techniques. 

In the following section, we outline the first-principles dynamical CPA presented in our paper II. After having 
established the basic formulation, it is desired how to calculate higher-order terms of individual harmonics in the 
dynamical part. In Sec. Ill, we calculate the higher-order terms using asymptotic approximation, and obtain the 
expressions for dynamical corrections. In Sec. IV, we present the results of numerical calculations for the densities 
of states (DOS) as the single-particle excitations in 3d transition metals. Calculations have been performed at high 
temperatures, where the present approach works best. We will demonstrate that the dynamical CPA can explain a 
systematic change of the DOS in 3d series from Sc to Cu. We also show that the screening effects on the exchange 
energy parameter are significant for the description of the excitation spectra in Mn, Fe, and Co. In the last section, 
we summarize the present work and discuss future problems to be solved. 



II. FIRST-PRINCIPLES TB-LMTO DYNAMICAL CPA 



We consider here a transition metal system with an atom per unit cell, and adopt the TB-LMTO Hamiltonian 
combined with a LDA-I-U Coulomb interactions as follows^. 

H = Ho + Hi, (1) 



iLa 



iLa 



iLjL'a 



(2) 



Hi — ^ ^'^^UonUm-tnUmi + ^ (Ul — -J)fliimnilm' — ^ JSilm ' SUm' ■ (3) 
i Til m>m' m>rn' 

Here e^^ in the noninteracting Hamiltonian Hq is an atomic level on site i and orbital L, /i is the chemical potential, 
tiLjL' is a transfer integral between orbitals iL and jL'. L = {l,m) denotes s, p, and d orbitals. alj^^ {aiLa) is the 
creation (annihilation) operator for an electron with orbital L and spin cr on site i, and JT-iLcr — a\j^^aiLu is a charge 
density operator for electrons with orbital L and spin cr on site i. 

In the Coulomb interaction term Hi, we take into account on-site interactions between d electrons {I — 2). 
(Ul) and J in Hi are the intra-orbital (inter-orbital) Coulomb and exchange interactions, respectively, fium (sum) 
with / = 2 is the charge (spin) density operator for d electrons on site i and orbital m. It should be noted that 
the atomic level in Hq is not identical with the LDA atomic level e^. The former is given by the latter as^i^ 
e° = — dE^j^j^/driiLa- Here UiLcr is the charge density at the ground state, -Elda ^ LDA functional to the 
intraatomic Coulomb interactions. 
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The free energy of the system T is written in the interaction representation as follows. 



Tr 



Texp 



(i/o(r) + ffi(r)) 



(4) 



Here /3 is the inverse temperature, T denotes the time-ordered product for operators. Ho{t) {Hi{t)) is the interaction 
representation of Hamiltonian Hq (Hi). 

We transform the interaction Hi (r) in the free energy into a one-body dynamical potential v making use of the 
Hubbard-Stratonovich transformationi2i2i. The transformation is a Gaussian formula for the Bose-type operator 



/ in^^"] ^ 



mm' '^m' 



(5) 



Here Amm' is a M x M matrix, and {a;,„} are auxiliary field variables. The above formula implies that the two-body 
interaction X)mm' bmAmm'bm' is transformed into a one-body interaction with a potential — X^m' '^^mm'Xm' coupled 
with the random fields {xm'}- 

After making use of the transformation at each time r, the free energy J- is written in the Matsubara frequency 
representation as follows. 



. N 21+1 

/ [n n ^^jrnSQn] CXp [ - /^i?^, C] ] , 



(6) 



(7) 



Here N is the number of sites, Cim(*w„) (^imQ(*'^Ti)) is the n-frequency component of an auxiliary field C,im{T) {£,ima{T)) 
being conjugate with inix{T) {rhiLaiT) — 2sii^cc{T)) for ^ = 2. Sp in the second term at the r.h.s. (right-hand-side) of 
Eq. ([7]) denotes a trace over site, orbital, frequency, and spin, g is the temperature Green function for noninteracting 
system _ffo- The matrices and are defined by A^m' = f/o<5mm' + (2C/i - J)(l -^mm'): ^mm' = Ji^-^mm') 

(a = X,y), and Bf„„, = UoSram' + - ^mTn')- 

The functional integrals in Eq. ([6]) are defined by 



. 21+1 



[ n 



N 



n fu^)2i+i n ^c™(o) 

-VI — 1 V ^ ' m — 1 



f^^det^^. 



21+ 



—d Cim(iw„) 



(8) 



Here (PQm{iujn) — 'iReCi„i(icj„)(iImCi„j(iw„). The dynamical one-body potential v at the r.h.s. of Eq. ([7]) is defined 
by 



iv)iLnajL'n'cr' = WjLcrcr' (jw„ - iuJn')SijS_ 



LL' 



(9) 



iL<j<T' {ii^^) = -77 i^mm'Cim' {i^u)Sl2Saa' - - ^ ^ -B.mm'Cjm'Q (iw^)<5;2 (^q ' 



(10) 



Ga {a = X, y, z) being the Pauli spin matrices. 

In the dynamical we introduce a site-diagonal coherent potential 

(S)iLncrjX'n'o-' = ^L(j{i^n)5ij5 LL' 5nn' ^crcr' , (11) 

into the potential part of the energy functional C]i and expand the correction w — E with respect to sites. 

E[i, C] = A^^(S) + E^[i^X^] + AE . (12) 
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Here the zero-th order term is a coherent part of the free energy which is defined by 



(13) 



The next term in Eq. (|12p is a sum of the single-site energies £'i[^i, Ci]- The dynamical CPA neglects the higher-order 
terms AE associated with inter-site correlations. 
The free energy per atom is finahy given by^SiSi 



(14) 



Here Jx = Jy = J± = [I — I / {21 + I)] J, Jz = Uo/{2l + 1) + J±, and we expressed the single-site term (the second 
term at the r.h.s. Eq. (IT4)) ) with use of an effective potential i?cff(C) projected onto a large static field variables 
— J2m Cmct(O). Morcover we have omitted the site indices for simplicity. 

The effective potential £'efi(^) consists of the static contribution i?st(^) and the dynamical correction term i?dyn(^)- 



The former is given as 



(15) 



(16) 



Here 6vLaiO) = wl(t(0) - TiLaii^^n), and ^ + ul^(O) is a static potential given by VLaiO) = [{Uq - 2Ui + 

J)nimi£,) + {2Ui — J)ni{^)]/2 — The electron number nL(^) for a given ^ is expressed by means of an impurity 

Green function as 



(17) 



and ni{^) = J2m '"■^(C)- The impurity Green function GLa{^, i^n) has to be determined self-consistently. The explicit 
expression will be given later (see Eq. ([5^ ). 

The coherent Green function Flo-(«w„) in Eq. (fT6|) is defined by 



FLaiiUJn) = [{i^n Hq S(iw„)) ^LaiLa- 



(18) 



Here (J?o)iL(TjL'(T is the one-electron Hamiltonian matrix for the noninteracting Hamiltonian Hq, and 

{T,{iu}n))iLajL'cr = ^La{i^^n)k]5LL' ■ 

The dynamical potential i?dyn(^) in Eq. ([15]) has been obtained within the harmonic approximation^Si^i^iSS. It is 
based on an expansion of i?dyn(C) with respect to the frequency mode of the dynamical potential VLaa'ii^v), where 
LJi, = 2vn / 13. The harmonic approximation is the neglect of the mode-mode coupling terms in the expansion. We 
have then (see Eq. (55) in our paper II) 



(19) 



Here the determinant Di, is a contribution from a dynamical potential vlo-o-' {'i-'-^u) with frequency uj^, and the upper 
bar denotes a Gaussian average with respect to the dynamical charge and exchange field variables, C,„(iw„) and 
Cma(iw„) (a = x,y,z). 

The determinant D^, is expressed as2i 



'21+1 



n D,.{k,m) 



(20) 
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D^{k,m) 



1 1 

a-„+kii^,rn) 1 

akiv,m) 







1 

1 1 

^k{v,m) 1 1 



(21) 



Note that 1 in the above determinant denotes the 2x2 unit matrix, a„(i^, m) is a 2 x 2 matrix defined by 



VLaa" {■li^u)9La"a"' {i^^n " iu)^)vLa"' a"" {-i^^)gLcr"" a' {i^n) , 



a (J a 



(22) 



(23) 



Here VLaa'ii^v) is defined by Eq. (|10p. gLaa'ii^n) is the impurity Green function in the static approximation, 
(Fl(iw„))o-o-' = FL„{iUn)5aa' , and 5wo is defined by {5va)„„, = vl„„,{0) ~ Y.L„{iUn)5aa' ■ 
The determinant D^{k,m) is expanded with respect to the dynamical potential as follows. 



D^{k,m) = I + D^J\k,m) + D';^\k,m) ^ , 



(24) 



(fc, to) = ^ -Wai (J^, to)w7i {-iy,m)---Va„ {v, to)u^„ (-i^, m)i)|^'^j {v, k, to) 



ai7i ■ ■ -071771 

Here the subscripts and 7, take 4 values 0, x, y, and z, and 



(i/,to) = -^i^ A,„„/C„/(iWi.)(5;2 



(25) 



(26) 



Va{i^,m) = -i^S"„„,^„i.a(iw^)(5/2 , (a = x,y,z) . 



(27) 



Note that the subscript {07} of i)|^^^j(i', A;, to) in Eq. p5|) denotes a set of (ai7i, • • • , a„7„). The frequency dependent 
factors l)|^^^j(j/, fc, to) consist of a linear combination of 2n products of the static Green functions. Their first few 

terms are given in Appendix A of our paper H-^ . Approximate expressions for higher-order terms will be given in the 
next section. 

Substituting Eq. (|24|) into Eq. (|20|) and taking the Gaussian average, we reach 



(«) 



n—1 



and 



my 



21+1 u-l 

E E Enn 

^l{k,rn)—n {aj{k,m)} P m—lk—0 



l{k,m) 

(n 

i=i 



(28) 



(29) 



Here each element of {l{k, m)} (fc = 0, • ■ • , — 1; to = 1, • • • ,21 + 1) has a value of zero or positive integer. aj{k, m) 
takes one of 4 cases 0, x, y, and z. j denotes the j-th member of the (fc, to) block with l{k, m) elements. P denotes a 
permutation of a set {(j, fc, to)}; P{(j, fc, to)} = {(jp, fcp, TOp)}. ctp-i means a rearrangement of {aj{k, m)} according 
to the inverse permutation P~^. The coefficient C^^, in Eq. (|29|) is a Coulomb interaction defined by 



(a = 0) 

(a = a;,y,2:) 



(30) 
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The coherent potential can be determined by the stationary condition SFcpa/STi = 0. This yields the CPA equation 
as^ 

{GLait^iO,,)) = FLai^^n) . (31) 

Here ( ) at the l.h.s. (left-hand-side) is a classical average taken with respect to the effective potential Ecs{C)- The 
impurity Green function is given by 

EE- 



GLa{£,, ii^l) = 9Laa{ii^l) + ^ ^ ■ (32) 

n f 

Note that the first term at the r.h.s. (right-hand-side) is the impurity Green function in the static approx- 
imation, which is given by Eq. (j23l) . The second term is the dynamical corrections, and KLa-{i^i) = 1 ^ 

Solving the CPA equation ([5T|) self-consistently, we obtain the effective medium. The electron number on each 
orbital L is then calculated from 

Tier 

The chemical potential /i is determined from the condition Ue = Here rie denotes the conduction electron 

number per atom. The magnetic moment is given by 



ncr 

In particular, the / = 2 components of magnetic moment are expressed as 

(m,) = (0 . (35) 
The above relation implies that the effective potential i?off(^) is a potential energy for a local magnetic moment ^. 



III. HIGHER-ORDER DYNAMICAL CORRECTIONS IN ASYMPTOTIC APPROXIMATION 

In our previous paper H^-'^, we took into account the dynamical corrections up to the second order in Eqs. (j28p and 
([32|) . We will obtain higher-order terms in this section using an asymptotic approximation. 

We note that the coupling constants B^^, = B^^, = J(l — 6mm') are much smaller than Amm' and -B^„/ because 
Uq and Ui ^ J. The latter condition is not necessarily satisfied for Sc and Ti. But J in these elements are small as 
compared with the d band width. Thus we neglect the transverse potentials, Vx{v, m) and Vy{iy, m) in the higher-order 
dynamical corrections. This approximation implies that a„(i^, m)^— cr = 0. The determinant D^{k,m) in Eq. (j20l) is 
then written by the products of the single-spin components as 

D^{k,m) = D^f{k,m)D^l{k,m). (36) 

Here Dy„{k,m) is defined by Eq. (PT|) in which the 2x2 unit matrices have been replaced by 1 (i.e., 1x1 unit 
matrices), and the 2x2 matrices a„{v,m) have been replaced by the 1x1 matrices an(y,m)aa- The latter is now 
given by 

a„(i^, to)ctct = '^Va{iy,m)v^{-iy,m)ha-faena{i',m), (37) 

a, 7 



e„^(j/, m) = gLa{n - v)gi^„{n). (38) 
Here ha-^^ = ^a-y + a(\ — Ja^), and we used a notation ghain) — gLaa-ii^n) for simplicity. 
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With use of the Laplace expansion, the determinant D^„{k,'m) can be written as 

D^a{k, to) = Dia{v, k, m)Dia{v, k, to) - akiv, 'm)aab2a{v, k, m)D2a{i^, k, to) 
Here Dna(y, k, m) and Dna(y, fc, m) are defined by 

1 1 



1 

«(n+l)i/+fc(l', "T,)cr<T 1 1 



(39) 



(40) 



DnA^, k,m) 



1 







a-niy+k(y,m)aa 1 1 

a-(n+2)i/+fc('^,TO)o-o- f f 







As we have shown in Appendix A in our paper I, Dna{i^, k, to) and I?„o-(z^, k, to) are expanded as foUows. 
Dna{v,k,m) = 



oo oo 



l + ^(-l)*^ X! ■'■ X! ai^v+k{v,m)ac,ai2v+k{v,m)c,a---ahv+k{v,m)c,cr, 

i=l li=nl2=li+2 (i=;,_i+2 



(41) 



(42) 



DnA^, k,m) 



OO oo 



i=l h=nl2=li+2 ii=ii_i+2 



Substituting Eq. ([37]) into Eq. (|42]), we obtain 
Dna{v.,k,m) = 

oo 0,^ 

X! X! '"ali'^'^)'"'fli~'^^^^)■■■'^a^i'^^'^)^"1^i~'^^^)^nU{o^l},'^,k,m) , 
i—1 Q:i7i ■ ■ ■a.i'^i 



(43) 



(44) 



-Di*2({a7}, fc, "i) = (-1)* /lai^K^ • • • K,~f,aA^^l{v, k, m) , 



(45) 



OO OO 



A^^Xi^^k.m) ^ ••• ^ e;ii.+fec,(z^, m)e;2^+fc<:r(i^,TO)---e;,i,+fc£r(i^,TO). 



In the same way, Dnai^, fc, to) is expressed by Eq. (|44p in which Ana{v, k, m) has been replaced by 



oo oo oo 



{ly, m)e^i^i,+ka{v, to) • • • e-i^u+ka{v, m). 

lx=nl2=li+2 ii=ii_i+2 



(46) 



(47) 



The quantities An^iv, k, to) and Anl^v, k, ni) contain the i-fold summations. In order to reduce these summations, 
we make use of an asymptotic approximation. The approximation is based on a high-frequency behavior of gLa{n) as 



9La{n) 



O 



(48) 
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Then the product of gLa{n — v) and gLain) in ena{i^, fn) is written by their difference as 

ena{i^, in) q^{gLa{n -f) - gLa{n)) , (49) 

where = P/2TTvi. 

Substituting Eq. into Eq. successively, we find 

A^^l{v,k,m) ^ — qlgLa{{n-l)v + k)gLa{nv + k)---gLa{{n + i-2)v + k). (50) 
I ! 

In the same way, we have 

A\ll{u, k,m) ^ —,ql gLai-nv + k)gLa{-{n + l)i^ + fc) • • -gLai-in + i - + k) . (51) 
Substituting Dna{>^ik,m) with Eq. (jSOp and Z)„o-(j^, fc, m) with Eq. (|5ip into Eq. p9p. we obtain 

CO 0,2 

D^.a(fc,TO) = ^ ^ Uai(i^,rn)u7i(-i^,m) • • •■!;a,(j^,m)?;,^,(-j/,m)L)[2-y}<^(f^,fc, w) , (52) 

1=1 ai7i---ai7i 

^g7>.(^>fc,"^) = AW({a7})|fi?(')(^^,fc,m) . (53) 



Here 



i-i 



B^j\v,k,m) = [\{~gLa{j^ + k) 



+ ES— v[ n 5L.(j^^ + fc)][l + ^5L.(*;^ + fc)] , (55) 

and D^^^^y^{i', k, m) = 1. n;({a7}) in Eq. ()54p is the number of {ai7i} pairs such that = 7^ among the I pairs. 
Substituting Eq. ^ into Eq. (031), we obtain i)|"'^j(i^, fc, to) in the asymptotic approximation. 



i=0 



Here we wrote the subscript at the r.h.s. explicitly to avoid confusion. Note that the values of ai and 7^ are limited 
to or z in the present approximation. When there is no orbital degeneracy, Eq. ([55)) reduces to the result of the 
zeroth asymptotic approximation in our paper 

In the actual applications we make use of the exact form up to a certain order of expansion in -D[,™\ and for 
higher order terms we adopt an approximate form (j56p . In this way, we can take into account dynamical corrections 
systematically starting from both sides, the weak interaction limit and the high-temperature one. 

IV. NUMERICAL RESULTS OF EXCITATION SPECTRA 

In the numerical calculations, we took into account the dynamical corrections up to the second order (n < 2) 
exactly, and the higher-order terms up to the fourth order within the asymptotic approximation. Summation with 
respect to u in Eqs. and was taken up to = 100 for n = 1 and 2, and up to = 2 for n = 3, 4. 
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When we solved the CPA equation (|3ip , we adopted a decoupling approximation to the thermal average of impurity 
Green function—, i.e.. 




GUqVW)Aei.),^ujn) . (57) 



Here we wrote the static exchange field ^ as (^2,^1) so that the decoupling approximation we made becomes clearer. 
The approximation is correct up to the second moment {i.e., (^^)) and allows us to describe the thermal spin 
fluctuations in a simple way. 

On the other hand, we adopted a diagonal approximation^Sdl to the coherent Green function at the r.h.s. of Eq. 

ID. 

f,,,„, = /. ^s."^;)-'' ^ . ,58, 

J «w„ - e - I]LCT(ia;„) - AeL 

Here p]f^{e) is the local density of states for the LDA band calculation, and Ae^ = (e^ — £l)'^'2- The approximation 
partly takes into account the effects of hybridization between different / blocks in the nonmagnetic state, but neglects 
the effects via spin polarization. 

The CPA equation with use of the decoupling approximation (|57p yields an approximate solution to the full CPA 
equation (|3ip . For the calculations of excitation spectra, one needs more accurate solution for the CPA self-consistent 
equation. We thus adopted the following average i-matrix approximationi^i (ATA) after we solved Eq. (I5T]) with 
the decoupling approximation. 

yATA/- , N y , {GLa{£.z,(i,ii^n)) - FLajiUJn) .^„^ 

(G-Lo- (f 2 , Cl , «W„) )i'Lcr («W„) 

Here the coherent potential in the decoupling approximation is used at the r.h.s., but the full average ( ) of the 
impurity Green function is taken. The ATA is a one-shot correction to the full CPA (PT|) . 

The coherent potential T,l^{z) on the real axis z — lu -\- id is then calculated by using the Pade numerical analytic 
continuation method--. Here 6 is an infinitesimal positive number. The densities of states (DOS) as the single-particle 
excitations, pl{'^) are calculated from the relation, 

Pi(w) = --lmFi,(z) . (60) 

TT 

We adopted the same lattice constants and structures as used by Andersen et. al^ in order to investigate a 
systematic change of excitations. For fee Fe, we used the lattice constant 6.928 a.u. being observed at 1440 K. 
The LDA calculations have been performed with use of the Barth-Hedin exchange-correlation potential to make the 
TB-LMTO Hamiltonian In the present work all the dynamical CPA calculations have been performed at 2000 K 
in the paramagnetic state. 

We adopted average Coulomb interaction parameters U obtained by Bandyopadhyav et. al^, and the average 
exchange interactions J obtained from the Hartree-Fock atomic calculation^. The intra-orbital Coulomb interaction 
Uq, inter-orbital Coulomb interaction Ui, and the exchange interaction energy parameter J were calculated from U 
and J as Ui) — U + 8J/5, Ui — U ~ 2J/5, and J = J, using the relation U^) — Ui + 2J. 

Calculated Coulomb interactions from Sc to Cu are plotted in Fig. 1 as a function of the valence electron number 
Ug. For Fe and Ni, we adopted the values used by Anisimov et. al^, which are also shown in the figure. Recent 
calculations suggest that the exchange interactions in the 3c? metals are reduced by about 30% as compared with their 
atomic values^. These values are also shown in Fig. 1 by dotted line. We will discuss the screening effects of J on 
the spectra using the values. 

Before we present the results of excitation spectra in 3d series, we briefly discuss the 4th-order dynamical effects. 
Figure 2 shows the d partial DOS for the paramagnetic bcc Fe on various levels of approximations. The static 
approximation with the decoupling scheme (j57p causes too strong thermal spin fluctuations with large exchange 
splitting, and yields the two-peak structure as shown by a thin curves in Fig. 2. The second order dynamical 
corrections suppress the thermal spin fluctuations and reduce the d band width as well as the dip at a; = —0.12 Ry. 
The 4th-order corrections enhance the two peaks. Finally the ATA correction (the best result in the present work) 
reduces the peaks and shifts them towards the Fermi level {i.e., the low energy side). 

We have reported in our recent paper— that the 4th-order dynamical corrections improve the magnetic properties 
of Fe and Ni. We obtained the Curie temperature Tq = 2070 K (1420 K) for Fe (Ni) in the static approximation. The 
second-order dynamical corrections lead to Tq = 2020 K (1260 K) for Fe (Ni). The 4th-order dynamical corrections 
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FIG. 1: Intraatomic Coulomb and exchange energy parameters as a function of the conduction electron number of 3d transition 
metals. These parameters are obtained from the band^ and atomic^ calculations. Averaged Coulomb interactions U: closed 
circles and thin curve, intraorbital Coulomb interactions Uo'- solid curve, interorbital Coulomb interactions Ui: dashed curve, 
Hartree-Fock exchange interactions J: closed squares and dot-dashed curve, and screened exchange interactions 0.7J: dotted 
curve. U, Uo, Ui, J recommended by Anisimov et. al.^^ are plotted by Q, +, x, and □, respectively, for Fe (rie — 8) and Ni 
(ne = 10). 
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FIG. 2: The d partial densities of states (DOS) of bcc Fe at 2000K in various approximations. The DOS in the static 
approximation with the decoupling scheme (|57p : thin solid curve, the DOS with the 2nd-order dynamical corrections and the 
decoupling (|57|l : dotted curve, the DOS with the 4th-order dynamical corrections and the decoupling (|57|) . and the DOS with 
the 4th-order dynamical corrections in the ATA. 
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FIG. 3: Densities of states (DOS) as the single-particle excitation spectra for fee Sc. DOS in the static approximation: thin 
solid curve, DOS with dynamical correction: solid curve, DOS in the local density approximation (LDA): dashed curve. The 
X-ray photoemission spectra (XPS)2^ and inverse photoemission spectra (BIS)^''' for hep Sc at room temperature are shown by 
dotted curves. Note that these data are arbitrarily renormalized in order to fit the calculated DOS. 



further reduce Tc to 1930 K for Fe and 620 K for Ni, respectively. The latter is in good agreement with the 
experimental value 630 K, while the former is still overestimated by a factor of 1.8. A large reduction of Tc in Ni due 
to the 4th-order dynamical corrections has been attributed to a reduction of the DOS at the Fermi level. Calculated 
effective Bohr magneton numbers in the 2nd-order dynamical CPA are 3.0 /xb and 1.2 /ib for Fe and Ni, respectively. 
The 4th-order dynamical corrections yield 3.0 ^b for Fe and 1.6 /itB for Ni, both of which are in good agreement with 
the experimental values 3.2 /iB (Fe) and 1.6 /iB (Ni). 

Among the 3c? transition metals, scandium has the weakest Coulomb interaction as shown in Fig. 1. Calculated 
DOS for fee Sc are presented in Fig. 3. The DOS below the Fermi level is close to the LDA DOS except some detailed 
structures. This is due to a small number of d electrons per orbital (~ 0.3 per d orbital) and rather weak Coulomb 
interactions. It should be noted that there is no correlation correction to the sp bands in the present approximation 
([58|) . so that spiky structure of the sp bands in the LDA calculations remains in the total DOS. The d DOS are 
smoothed by the scattering corrections of the self-energy. The corrections become larger near the top of d bands, so 
that the peak of t2g band at w = 0.275 Ry is much damped down, the d band becomes narrower, and the spectral 
weight shifts to the higher energy region. The difference between the dynamical and static DOS is rather small except 
high energy region (w > 0.2 Ry). We performed the numerical calculations with use of the screened value J — 0.029 
Ry, but we hardly found the change of DOS in shape. 

Calculated DOS qualitatively agree with the XPS and BIS data for hep So^iiS^ as shown in Fig. 3. (Note that the 
crystal structure of Sc is not the fee but the hep experimentally.) Here and in the foUowings the intensities of the 
experimental data are arbitrarily scaled to fit theoretical DOS. Rapid decrease of the XPS and BIS data indicates 
the cut-off due to Fermi distribution function. Moreover the deviations from the DOS in high energy region are due 
to secondary electrons, and outside the scope of the present theory. The high-energy peak in the calculated DOS 
around oj = 0.25 Ry deviates from the BIS peak at a; = 0.30 Ry. This is partly explained by the difference in crystal 
structure. In fact, the LDA calculations— indicate that the peak position in the hep Sc is higher than that of the fee 
one by 0.025 Ry. A small hump at about oj = —0.2 Ry in the XPS data does not appear in the present calculations. 
The discrepancy is not due to the hep crystal structure since there is no corresponding peak in the LDA DOS for the 
hep Soil. 

In the case of fee Ti, we obtained the DOS being similar to the fee Sc as shown in Fig. 4. Thermal excitations 
smooth the LDA DOS, damp the highest peak at w = 0.25 Ry, and transfer the spectral weight to the higher energy 
region (w ^ 0.30 Ry). The dynamical corrections are not so important in the case of Ti as understood by comparing 
the DOS with that in the static approximation. We find rather good agreement of the DOS with both the XPS and 
BIS experimental data for hep Ti. 

We have calculated the excitation spectra of vanadium for the bee structure as shown in Fig. 5. In this case the 
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FIG. 4: Calculated DOS for fee Ti. The notations are the same as in Fig. 3. The XPS^'* and BIS^ data for the hep Ti are 
obtained at room temperature. 
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FIG. 5: Calculated DOS for bee V. The notations are the same as in Fig. 3. The XPS^i and BIS2^ data for the bcc V are 
obtained at room temperature. 

crystal structure is identical with the experimental one. The main peaks and valleys in the LDA DOS are much 
weakened by local electron correlations, and the spectral weights move to the higher energy region. The d bands in 
the quasiparticle energy region < 0.2 Ry) shrink by about 10% as compared with the LDA one. The calculated 
DOS shows a good agreement with the XPS and BIS dataS^!^ in lineshape. Note that any artificial parameter is not 
introduced for comparison between the theory and experiment. The DOS in the static approximation yields an excess 
d band broadening due to thermal spin fluctuations. 

Calculated excitation spectra of the bcc Cr is similar to that in the bcc V as shown in Fig. 6. Because the valence- 
electron number of Cr is larger than that of V by one, the Fermi level shifts to the higher energy region. The t2g peak 
around oj = —0.15 Ry in the LDA DOS is weakened, and shifts toward the Fermi level. The position of the Cg peak 
at w = 0.1 Ry is not changed, but its weight is much decreased by electron correlations. The static approximation 
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broadens the Cg peaks excessively, thus does not explain the BIS data. Calculated peak around w — —0.1 Ry seems 
to be too small as compared with the XPS experiment. There is a possibility that the antiferromagnetic correlations 
enhance the peak in the present calculations, because the experimental data are taken below the Neel temperature. 

The Coulomb interactions of Mn are roughly twice as large as those in Sc, while their LDA DOS in the fee structure 
are similar to each other. In addition, electron number per atom increases from 3 to 7. Thus one expects more 
electron correlations in the case of Mn. Calculated DOS as well as experimental XPS-BIS data^'*'^^ are shown in Fig. 
7. We find that the central peak consisting of the t2g bands around to = —0.1 Ry in the LDA DOS changes to a 
valley due to electron correlations, so that the DOS shows a two-peak structure. The result indicates a formation of 
the Mott-Hubbard type bands due to strong on-site correlations as we suggested in our recent paper—. The same 
result was obtained by Birmann et. al.^^ by using the Hamiltonian without transverse spin fluctuations. Static 
approximation overestimates the splitting of the Mott-Hubbard bands. The dynamical effects suppress such a band 
broadening due to thermal spin fluctuations. When we adopt the screened value J = 0.043 Ry, the DOS around 
uj — —0.1 Ry increases and the two-peak structure becomes less clear. 

Experimentally, the bulk fee Mn is realized only in the narrow temperature range between 1352 K and 1407 K at 
high temperatures. There is no photoemission experiment in this temperature regime. However the XPS data for 
the 20 monolayer fee Mn on Cu3Au(100)^ are available at room temperature. The data seem to be explained by the 
dynamical CPA with partially screened J between 0.043 and 0.061 Ry. Theoretical results agree with the BIS dataiS^ 
in peak position. 

The line shape of the calculated DOS for bcc Fe is similar to the bcc Cr as shown in Fig. 8. The main peak of the 
LDA DOS near the Fermi level and the peak of t2g bands around oj = —0.15 Ry are much weakened due to electron 
correlations. The spectral weight moves to higher energy region. The static approximation overestimates the band 
width by 20%. The main peak at u = 0.04 Ry above the Fermi level is consistent with the BIS datai^ at 0.86Tc, 
Tc being the Curie temperature. (Note that a hump at a; = 0.1 Ry in the BIS data is the remnant of the e2g peak 
for the minority band and should disappear above Tc-) On the other hand, the DOS below the Fermi level does not 
well correspond to the XPS dataSi at l.OSTc. One needs more weight around oj = —0.1 Ry in order to explain the 
experimental data. 

The DOS calculated with use of the screened value J = 0.046 Ry considerably enhance the weight around w = —0.1 
Ry, thus the screening on J partly explains the broad peak at lj — —0.1 Ry in the XPS data. This feature does not 
change even if we adopt U value obtained by Bandyopadhyav et. al. The magnetic short-range order— may also 
explain the discrepancy between the XPS data and the present result based on the single-site approximation because 
the experimental data of the XPS are taken near Tq- 

We have also investigated the DOS for the fee Fe. The fee Fe is well-known to be a typical itinerant magnet showing 
the spin density waves with magnetic moment of about 1 /ie per atom at low temperatures^^—. But the fee Fe shows 
anomalous thermal expansion^®, so that the d band width at high temperatures is expected to become narrower than 
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FIG. 7: Calculated DOS for fee Mn. The notations are the same as in Fig. 3. The dot-dashed curve is the dynamical results 
for U = 0.192 Ry and the screened exchange energy parameter J = 0.043 Ry. The XPS datai^ for the fee thin-film Mn on 
Cu3Au(100) and the BIS data^^ for a-Mn at room temperature are plotted for comparison with the theory. 
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FIG. 8: Calculated DOS for bcc Fe. The notations are the same as in Fig. 3. The dot-dashed curve is the dynamical results 
for U — 0.169 Ry and the screened exchange energy parameter J — 0.046 Ry. The XPS data^ are measured at 1.03Tc, Tc 
being the Curie temperature. The BIS dataiS^ are measured at Q.86Tc. 



those at low temperatures by several percent. We may then expect stronger correlation effects. As shown in Fig. 9, 
we find that the DOS calculated with use of the unscreened value J = 0.066 Ry is similar to that of the fee Mn. The 
t2g central peak in the LDA DOS splits into the lower and upper Mott-Hubbard bands due to on-site correlations. 
The band splittings are smaller than those in the fee Mn, and the dip at cj = —0.15 Ry is weakened (see Fig. 7). The 
static approximation overestimates the band width by 20 %. 

We present in Fig. 10 the d partial DOS for Cg and t2g bands in order to examine the details of the band splitting. 
As seen from the figure, both the Cg and the t2g DOS show the two-peak structure. The energy difference between 
the upper and lower peaks is about 0.25 Ry, which is approximately equal to the intra-orbital Coulomb interaction 
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FIG. 9: Calculated DOS for fee Fe. The notations are the same as in Fig. 3. The dot-dashed curve is the dynamical results 
for U = 0.169 Ry and the screened exchange energy parameter J = 0.046 Ry. The XPS^* and BIS£2' data are obtained for the 
fee Fe on Cu(lOO) at room temperature. 
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FIG. 10: Partial d DOS of fee Fe for eg (dashed curve) and t2g (solid curve) electrons. The total d DOS is shown by thin solid 
curve. The partial DOS in the LDA are also shown by dotted curve (cg), dot-dashed curve (t2g), and thin dot-dashed curve 
(sp electrons), respectively. 



for fee Fe, Uq = 0.27 Ry. Moreover, we find a large scattering peak of — ImSLo-(w -I- iS) at w = —0.12 Ry. These 
behaviors verify that the two-peak DOS forms the Mott-Hubbard bands due to electron correlations. 

It is remarkable that the dynamical results are sensitive to the choice of U and J in the case of fee Fe. The two-peak 
structure almost disappears when we adopt U = 0.169 Ry and the screened value J = 0.046 Ry, as shown in Fig. 9, 
while it again appears when we adopt U = 0.219 Ry obtained by Bandhyopadhyav it et. al. and the screened value 
J — 0.046 Ry. The behaviors are understood from the following arguments. The Coulomb interactions suppress the 
doubly occupied states of electrons on an orbital in general. The Hund-rule coupling J also suppresses the doubly 
occupied states to reduce the energy; J tends to enhance U effectively. Thus the increment of C/ or J is favorable for 
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FIG. 11: Calculated DOS for fee Co. The notations are the same as in Fig. 3. The dot-dashed eurve is the dynamical results 
for U = 0.245 Ry and the screened exchange energy parameter J — 0.048 Ry. The XPS^ and BIS^^ data for the hep Co at 
room temperature are drawn by the dotted curves. 



the formation of the Mott-Hubbard bands. 

The peak near the Fermi level explains well the BIS data^^ for the fee Fe on Cu(lOO). The XPS data^^ for the fee 
Fe on Cu(lOO) are somewhat controversial. The peak around uj = —0.2 Ry is usually interpreted as a peak due to Cu 
substrate. If this peak originates in the Cu substrate by 100%, the data support the result for {U, J) = (0.169,0.046) 
Ry. But if the peak is interpreted as a superposition of both the fee Fe and the Cu substrate spectra, we expect 
a two-peak structure of fee Fe, and the XPS data are consistent with the results for (J7, J) = (0.169,0.066) and 
([/, J) = (0.219,0.046) Ry. Resolving the problem is left for future investigations. 

The DOS of fee Co in the paramagnetic state does not show the Mott-Hubbard type structure any more as shown 
in Fig. 11. The peak of the t2g bands at a; = —0.15 Ry becomes almost flat, and the peak of the t2g bands on the 
Fermi level is much weakened. We find a weak hump around oj = —0.5 Ry suggesting the '6 eV satellite. When we 
adopt the screened value J = 0.048 Ry, the peak at w = —0.15 Ry partially remains as found in Fig. 11, and the 
hump around lo = —0.5 Ry almost disappear. 

We present in Fig. 12 the Cg and t2g partial DOS to clarify the disappearance of the two-peak structure in the total 
DOS. Calculated t2g partial DOS does not show the band splitting any more. The Cg DOS still show the two-peak 
structure. But the dip between the peaks is shallower than that of the fee Fe (see Fig. 10) and the energy difference 
between the peaks (0.18 Ry) is smaller than the intra-orbital Coulomb interaction energy Uq = 0.36 Ry. In the case 
of Co, the number of hole states is reduced as compared with that of fee Fe. This reduces the charge fluctuations 
and thus the magnitude of the self-energy by a factor of two. The reduction yields the disappearance of the two-peak 
structure in the total DOS. 

There is no experimental data on the fee Co above the Curie temperature. Calculated DOS does not show a good 
agreement with the XPS^^ and BIS data^^ of hep Co at room temperature. The UPS (Ultraviolet Photoemission 
Spectroscopy) data for hep Co by Heimann et. al.^^ show the existence of the 6 eV satellite in agreement with the 
present result for the unscreened J, while the other dat a^'^i^" do not. 

Single-particle excitations of Ni have been investigated extensively in both theory^iiStiS and experiment^"—. 
Present result of fee Ni shows a single-peak structure as shown in Fig. 13. Moreover the correlations increase the 
spectral weight around uj = —0.45 Ry, and creates a small hump corresponding to the 6 eV satellite due to two- 
hole excitations. The dynamical effects suppress the band broadening of the static approximation by 20 %. These 
behaviors do not change even if we adopt the screened value J = 0.046 Ry. The screened J enhances the main peak 
around uj = —0.05 Ry and creates a hump at a; = —0.15 Ry in the DOS. Though the calculated DOS is consistent 
with the XPS^ and BIS data^-, the band width seems to be somewhat larger than that of the XPS data. 

We present finally in Fig. 14 the DOS for Cu. Electron correlations via a small number of d holes (ss 0.36 per atom) 
move the spectral weight of the LDA DOS to the lower energy region. The peak of the d bands shifts toward lower 
energy by 0.05 Ry. It is also remarkable that a broad hump appears at w = —0.7 Ry due to two-hole excitations. We 



17 




Energy (Ry) 

FIG. 12: Partial d DOS of fee Co for eg (dashed eurve) and t^g (solid eurve) eleetrons. The total d DOS is shown by thin solid 
eurves. The partial DOS in the LDA are also shown by dotted curve (eg), dot-dashed curve (t2g), and thin dot-dashed curve 
(sp electrons). 
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FIG. 13: Calculated DOS for fee Ni. The notations are the same as in Fig. 3. The XPS^ and BIS^'^ data for the fee Ni are 
obtained at room temperature. 

find a good agreement between the dynamical CPA theory and the experimen t^^'^^ for this system. 

V. SUMMARY 

In the present paper, we have obtained approximate expression of the higher-order dynamical corrections to the 
effective potential in the first principles dynamical CPA, making use of an asymptotic approximation. The approx- 
imation becomes exact in the high frequency limit, and much reduces the multiple summations with respect to the 
Matsubara frequency at each order of expansion in the dynamical corrections. 
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FIG. 14: Calculated DOS for fee Cu. The notations are the same as in Fig. 3. The XPS ■ and BlS^a data for the fee Cu are 
obtained at room temperature. 

Within the 4-th order dynamical corrections, we have investigated systematic change of the DOS at high tempera- 
tures in 3d series from Sc to Cu. Thermal spin fluctuations in the static approximation smooth the LDA DOS at high 
temperatures, especially reduce their main peaks and broaden the d band width. Dynamical effects reduce the band 
broadening, and move the spectral weight to higher energy region. These effects explain in many cases the lineshape 
of the XPS and BIS experimental data from Sc to Cu quantitatively or semiquantitatively. We found the formation 
of the Mott-Hubbard type bands due to electron correlations in the case of fee Mn and fee Fe, and also found that 
the dynamical effects can create a small hump corresponding to '6 eV satellite in Co, Ni, and Cu. 

We investigated the effects of the screened exchange energy parameters using the reduced values of J by 30 %. The 
screening of J is significant for the DOS in Mn, Fe, and Co. The reduction of J tends to weaken the Mott-Hubbard 
type bands in fee Mn, and even could destroy the two-peak structure in the case of fee Fe. It also develops the 
central peak around uj — —0.1 Ry in bcc Fe and fee Co. Some of these results explain better the XPS data. But we 
have to calculate the other physical quantities with use of the same scheme, and have to examine in more details the 
consistency among them in order to conclude the validity of the screened values of ,/. The magnetic short-range order 
should also be important for understanding the experimental data of magnetic transition metals for more detailed 
discussions. 

Present theory explains a systematic change of the spectra in 3d series at high temperatures. At lower temperatures, 
the higher-order dynamical corrections should be more important. Improvements of the theory in the low temperature 
region are left for future investigations. 



Numerical calculations have been partly carried out with use of the Hitachi SRI 1000 in the Supercomputer Center, 
Institute of Solid State Physics, University of Tokyo. 
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